The relation between Feynman cycles and off-diagonal long-range order 
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The usual order parameter for the Bose- Einstein condensation involves the off-diagonal correlation 
function of Penrose and Onsager, but an alternative is Feynman's notion of infinite cycles. We 
present a formula that relates both order parameters. We discuss its validity with the help of rigorous 
results and heuristic arguments. The conclusion is that infinite cycles do not always represent the 
Bose condensate. 
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Introduction 

In 1953 Feynman suggested an order parameter for 
the Bose-Einstein condensation of interacting systems 
[1|. It involves the lengths of space-time trajectories in 
the Feynman-Kac representation. Three years later Pen- 
rose and Onsager introduced the concept of off-diagonal 
long-range order |2|. The latter is now accepted as the 
correct order parameter. Yet Feynman cycles remain ac- 
tual because they are easier to deal with, and because 
they provide a useful microscopic heuristic [j, H, |5| . An- 
other order parameter that involves space-time trajecto- 
ries winding in spatial directions — rather than in the 
imaginary time direction — is closely related to Feyn- 
man cycles; it has been used extensively in investigations 
of the supersohd phase, see 0, 01 a-nd references therein. 

The theory of the ideal gas has led to the expression 
"Bose condensate" to designate the particles that occupy 
the zero Fourier mode. The concept extends to inter- 
acting systems via the off-diagonal correlation function 
(T{x,y). The density of the Bose condensate is then 






lim 

\x-y\^ 



(T{x,y). 



Feynman's approach introduces the density of infinite cy- 
cles, q{oo), to be precisely defined below. The natural 
question is whether these densities are identical. 

The links between Feynman cycles and off-diagonal 
long-range order were explored by Siito, who proved that 
infinite cycles occur in the ideal gas below the critical 
temperature [S|, and that they do not occur above it i^. 
These results were later extended to the mean-field Bose 
gas [13, El Articles [1111111111111111113 all assume, 
or conjecture, that Pq = q{oo) in any Bose system with 
reasonable interactions. We will see, however, that this 
cannot be true in general. 

The purpose of this letter is to present an explicit re- 
lation between the off-diagonal correlation function and 
the densities of Feynman cycles. Namely, 

(^(.x,y)^ ^ Cn{x~y)Q{n)+Cac{x-y)Q{oo). (1) 

n > 1 

Here, Cn{x) are coefficients, and Q{n) denotes the density 



of particles in cycles of length n. We expect that 
^ Cn{x) ^ 1, and lim c„(x) — Caa{x), 

n — *oo 

for any fixed x. The latter statement is not obvious. 
Notice that the analogous statement for Q{n) is false in 
general: We always have lim„ Q{n) — 0, but it is possible 
that ^(oo) 7^ 0. In addition, wc should have 



\x\ 



lim Cn{x) — 0, 



(2) 



for any fixed n. But Coo{x) may converge to a strictly 
positive constant c. The formula and these properties 
can be established in the case of the ideal gas in the 
canonical ensemble [l^l- It is found that 



M 



^ g-xV4«/3 
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and 



Q{n) 
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q{oo) = max(0,p- pc)- 



Here, fi is the chemical potential (that depends on the 
density p), and pc is the critical density of the ideal gas. 
Notice that c = 1. 

The interacting gas constitutes a formidable challenge 
to theoretical physicists. The sole mathematical proof 
about the occurrence of Bose-Einstein condensation deals 
with the hard-core lattice model at the symmetry point 
|l3l |lj| . A weakly interacting system is widely expected 
to display a Bose-Einstein condensation at low tempera- 
ture. However, the question of whether interactions in- 
crease or decrease the critical temperature is still cur- 
rently debated; the discussion in |j] illustrates this point 
by quoting several references that draw contradictory 
conclusions. In recent years progress has been made in 
understanding Bogolubov's theory [l^- Remarkable re- 
sults dealing with the ground state of low density systems 
are reviewed in |lfil |. 

Two results can be rigorously derived that add cre- 
dence to the formula (^. First, for /i < 0, for any (3, and 
for any repulsive interactions, we have 



lim (j{x,y) 



0, and ^(oo) 



0. 



Second, a rigorous cluster expansion allows to establish 
Eq. |(2l in the regime where /it < 0, and where /S is small 
with respect to the interactions. These results are proved 
in [13. 

In this letter we provide a heuristic discussion of inter- 
acting systems, and we will conclude that 

* Po ~ ^(oo) when interactions are weak; 

• Pq = 0, q{oo) > in a crystal at sufficiently low 
temperature. 

It is not clear whether a regime of parameters exists 
where the constant c differs from and 1. 

Throughout this letter, we denote finite volume expres- 
sions in plain characters, and infinite volume expressions 
in bold characters. 



Feynman-Kac expressions 

In this section we introduce the relevant mathematical 
expressions in the Feynman-Kac representation of Bosc 
systems. We consider a system of identical bosonic parti- 
cles interacting with a pair potential. We study this sys- 
tem in the grand-canonical ensemble using the Feynman- 
Kac representation. The present discussion is necessarily 
brief, and we refer to 12j| for a more extended mathe- 
matical description. 

The grand-canonical partition function in a domain D 
(a d-dimensional cube of size L, volume V — L"^, with 
periodic boundary conditions) can be expressed using the 
Feynman-Kac representation [17|. With /3 the inverse 
temperature and /i the chemical potential, it is given by 
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(3) 



1 ^ i<j ^ k 



This expression is illustrated in Fig.^ In words, the par- 
tition function involves a sum over the number k of space- 
time cycles. For each cycle j, we sum over the winding 
number uj; we integrate over the initial position Xj] we 
integrate over a Brownian trajectory loj that starts and 
ends at Xj, winding hj times around the time direction. 
Here, dW^y{uj) is the Wiener measure for a trajectory 
UJ : [0, (3] ^' D such that a;(0) — x and w(/3) — y. Given a 
trajectory to with winding number n, we denote by U{uj) 
the interactions between its legs. That is, 



U{lu) = 



1 



^ i<j ^ n- 



U{uj{il3 + s) - uj{jl3 + s))ds. 



and n') have interactions 

Uiu;,u;')^ J2 -^ / U{u;{if3 + s)-uj'{jf3 + s))ds. 

^ i ^ n'-l 

In the two equations above, U{x) denotes the pair in- 
teraction potential between two particles separated by a 
distance |a;|. We suppose that U{x) is nonnegative and 
spherically symmetric. 




FIG. 1: The Feynman-Kac representation of the partition 
function for a gas of bosons. The horizontal plane represents 
the d spatial dimensions, and the vertical axis is the imagi- 
nary time dimension. The picture shows a situation with five 
particles and two cycles, with winding numbers 4 and 1. 

The density of particles in cycles of length n is given 
by the expectation of the "observable" -^ X]j=i "■'^"j,"- 
Some computations yield the following finite volume ex- 
pression 



g{n) 



— gPp-n 
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(4) 



Here, Z{uj) is the grand-canonical partition function for 
a system where the trajectory w is present. It is given 
by the same formula as Eq. ^, but with the additional 
term 



k 
He 
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pU(i. 



For repulsive interactions we always have Z[uj) ^ Z , with 
equality in the ideal gas. When the volume is finite we 
have 



n > 1 



where p is the density of the system, which depends on 
the chemical potential. Let Q{n) denote the thermody- 
namic limit of Q{n). Limits and infinite sums do not 
always commute because of a possible "leak to infinity" . 
Fatou's lemma garantees that 



Different trajectories uj and lo' (with winding numbers n 



E ^(") ^ p- 



n> 1 



This suggests to define the density of particles in infinite 
cycles by 



£.(00) = P- ^ Q{n). 



(5) 



n> 1 



The ofF-diagonal correlation function of Penrose and 
Onsager also has an expression in the Feynman-Kac rep- 
resentation. It involves an open trajectory that starts at 
X and ends at y, possibly winding several times around 
the time direction. Precisely, we introduce 



<7(x,y) 



E 

n > 1 



e^^" I dW:l\Lo) e~^^^^^ ^. (6) 



True to our convention we denote the thermodynamic 
limit by cr{x, y). 

We immediately obtain a finite volume relation be- 
tween g{n) and a{x,y). It is enough to consider a{0,x) 
because of translation invariance (we use periodic bound- 
ary conditions). From Q) and ©, we have 



a{0,x) 



E 

n> 1 



Cnix) g{n), 



(7) 



where the coefficients c„(x) are given by 



Cn{x) 



rnP 



dW^oT(^)e- 



-I3U(lo) 



Z{uo) 



CMe-^^^-^^ 



(8) 



We need to understand the thermodynamic limit of l|7|l. 
"Leaks to infinity" occur, and we cannot simply replace 
finite volume terms by infinite volume ones. The ideal 
Bose gas yields exact expressions and therefore consti- 
tutes a useful first step. As F — > cxd, Eq. Q becomes 

El 



(t{0,x) 



-x^/inl5 



n> 1 



Q{n) + q{co). 



(9) 



The latter term represents the interactions between uj 
and the configuration of trajectories lo. The environ- 
ment typically consists of trajectories spread all over the 
domain, with no significant fluctuation of density. When 
interactions are weak, the typical trajectory from to 2; 
does not depend much on the presence of the random en- 
vironment. For n large we have U{uj,uj_) « un for almost 
all Lo and w, where m is a positive number. It is known 
that (weakly) self-avoiding random walks satisfy a local 
central limit theorem in high dimensions [l^illSl) so that, 
as n — > CO, 



C^jx) 
C^(0) 



exp 



x^ \ 
Aan/sJ' 



(10) 



The constant a is greater than 1 and it reflects the (weak) 
self- avoidance of to. The situation is then similar to the 
ideal gas, and we expect the formula |^ to hold with 
Cn{x) given by l(T?Hl . and Coo{x) = 1 for all x. 



Feynman cycles in crystals 

Bosons with strong interactions are likely to undergo 
a usual condensation into a solid at low temperature. 
We consider a crystalline phase where particles display 
(diagonal) long-range order, and where there is no off- 
diagonal long-range order. (We do not consider a super- 
solid phase here.) We are going to argue that infinite 
cycles are present at very low temperature, in dimension 
three or more. The typical Feynman-Kac representation 




Weakly interacting systems 

We seek to understand the general behavior of the co- 
efficients c„ and Coo m a gas of bosons with weak inter- 
actions. We can represent the numerator and denomina- 
tor of (|HJl as the connectivity of a weakly self- avoiding 
random walk in a random environment. The random 
environment is due to Z{ll!), which involves sums and in- 
tegrals over configurations lo of closed trajectories. We 
obtain 



diy;/(c.)e-''^('' 



Z{^) 



= M^)C^{x), 



with V a probability measure (Jdi^(a;) = 1), and Cn{x) 
the connectivity of the random walk between positions 
and x: 



C^{x)^ jdW^i{i^)i 



-^puiuj) -t3U{u,uy) 



FIG. 2: Schematic typical space-time configuration for a crys- 
talline phase. Trajectories are roughly vertical, with occa- 
sional tunneling to a neighbor. 

for a quantum system in a crystalline phase is depicted 
in Fig. 12 Particles are represented by rather straight 
trajectories that are given by Brownian motions in the 
presence of an effective harmonic trap, due to the neigh- 
boring particles. Occasionally, two neighboring particles 
may exchange their positions through quantum tunnel- 
ing. Thus Feynman cycles are associated with the follow- 
ing simple dynamical system. Particles occupy the sites 
of a lattice, and to a each bond is associated a Poisson 
process that exchanges the corresponding particles. The 
rate is proportional to /? — the actual constant depends 
on the potential barrier for switching two particles, and 
on the entropy of Brownian paths that do it. Any site 
belongs to a Feynman cycle of arbitrary length. 



The cycle that contains a given site can also be viewed 
as a continuous-time discrete random walk. From time 
to f3, jumps on neighboring sites occur with a rate that 
is essentially independent of (3. If the position at time /3 
is identical to the original one, the cycle closes and has 
winding number 1. The position is otherwise again given 
by a random walk, but with the following two restric- 
tions: (a) No jump can occur if the resulting space-time 
picture has more than one particle at a given site; (b) if 
the random walk meets the other end of a bond, it has to 
jump in the opposite direction. (This process appears in 
the study [20j| of the Heisenberg ferromagnet.) It is illus- 
trated in Fig. 131 The first restriction is a self-avoidance 
condition; the second restriction amounts to an effective 
attraction. These should not be too important in dimen- 
sions higher or equal to three, where random walks are 
transient. 




FIG. 3: The Feynman cycle in a solid resembles a discrete 
random walk. The thick bond is crossed in both directions. 

This suggests that the typical Feynman trajectories for 
bosons in a crystal (in 3D) is infinite if /3 is large enough, 
and lim^^co q{oo) — p. On the other hand, the off- 
diagonal correlation function shows exponential decay in 
a crystal, and Pq = 0. 

Finally, let us understand the behavior of the coeffi- 
cients c„ (x) , Coo (a;) in a- crystalline phase, so as to con- 
vince ourselves that Eq. Q and the subsequent prop- 
erties are still valid. For simplicity, we assume a cubic 
structure, but the argument is more general. We inves- 
tigate the typical space-time configurations, which con- 
tribute the most to the numerator and denominator of 
Eq. © . It helps to understand first the situation without 
tunneling, and then to take it into account. The denomi- 
nator typically involves crystalline space-time configura- 
tions as in Fig. |21 The tunneling must be such that the 
particle located at the origin belongs to a trajectory with 
winding number n. 

The numerator of ^ is more complicated. One can 
identify two alternative typical behavior, that are de- 
picted in Fig.^(a) and (b). Both alternatives yield the 
same conclusion, namely that c„(a;) « e""'^' provided n 
is large (or infinite). This result should be contrasted 
with the weakly interacting case, where c„(a;) « 1 if 




(a) 




FIG. 4: The two possibilities for the typical open trajectory 
from to X. For clarity, occasional tunneling between verti- 
cal trajectories is not shown, unlike in Fig. |5| (a) An extra 
trajectory is superimposed on the crystalline structure, (b) 
The open trajectory is part of the crystalline structure. 



The first alternative involves an open cycle that is su- 
perimposed to the background crystalline configuration, 
which is left essentially undisturbed; see Fig.0](a). Let 
m be the winding number in absence of tunneling. The 
contribution of all trajectories is approximately 



(47r/3m) 



-d/2 Bfi'm -x^/iml3 



(11) 



Here, p* = /i — {U{uj,oj_)) is the effective chemical poten- 
tial. It is negative since the crystalline phase is stable. 
The maximal contribution comes from m « „,,, , , and 

the dominant term in Ullfl is e~^~'^*l^l . For n ^ \x\ 
one can argue that tunneling contributes equally to the 
numerator and the denominator of ((SJ, and therefore 



.(^) 



-fj,'\x\ 



The second alternative is depicted in Fig. 01(b). The 
open trajectory is part of the crystalline structure. Since 
the distance between particles is p^^'^, the open trajec- 
tory involves |x|/p'^/'' "one-directional tunnelings" . Each 
contributes a small factor to the ratio in IjHIl. Since tun- 
neling affects equally the numerator and the denomina- 
tor, we again find that c„(x) decays exponentially with 
\x\, uniformly in large n. 

As a consequence, Eq. Q) predicts that cr{0, x) ^ as 
|a:| —f oo, while q{oo) ^ 0. 

Conclusion 

We have proposed an exact relation, Eq. (^, between 
Feynman cycles and the off-diagonal correlation function. 
This relation involves the coefficients c„ which can be 
computed in the case of the ideal gas. It is also backed by 



a few rigorous results valid for interacting systems. Our 
study suggests that both order parameters are identical 
in the weakly interacting Bose gas, but not in strongly 
interacting systems in a crystalline phase. The density of 
particles in infinite cycles, ^(oo), may be strictly positive 
while ^Q = 0. An interesting open question is whether 
^(oo) > Qq > for a certain choice of parameters. 

I am grateful to the referees, whose comments allowed 
me to improve the clarity of this letter. 
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